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TECHNICAL PUBLICATION 


INTEGRATED PULSE DETONATION PROPULSION 
AND MAGNETOHYDRODYNAMIC POWER 

1. INTRODUCTION 


Contemporary interest in pulse detonation engines (PDE’s) for various aerospace propulsion 
applications arises from the increased performance potential associated with unsteady flow effects and 
from anticipated cost reductions associated with a highly simplified design configuration. Unfortunately, 
the unsteady operational attributes of pulse detonation engines complicate their analysis when including 
nonideal effects, and it is difficult to make an accurate assessment of representative performance and cost 
benefits. Furthermore, design and optimization strategies are at an early stage of development, and there 
are critical operational issues, such as overall engine power requirements, that are in need of thoughtful 
attention. These uncertainties must be adequately addressed before the PDE will ever be realized in practice. 

In general, PDE’s can either collect air from the surrounding atmosphere to serve as the oxidizer or 
they can carry their oxidizer on board. In the former case, PDE’s represent an attractive low-cost alternative 
to turbojets based on their potential for surpassing intrinsic limits on scalability, operating range, efficiency, 
and manufacturing costs. In the latter case, the device is referred to as a pulse detonation rocket engine 
(PDRE). The most common oxidizer for a PDRE is oxygen (0 2 ), which yields a combustible mixture that 
is easier to detonate than mixtures using air. Hydrogen (H 2 ) is probably the optimum fuel for PDRE’s 
because of its ease of detonation and its high performance; however, hydrocarbon fuels are also of 
considerable practical interest. 

PDE’s produce thrust by intermittently detonating a charge of fuel and oxidizer placed inside a tube 
with one end closed and the other end open. The resulting propagation speed of the combustion detonation 
wave is extremely high, implying orders of magnitude increase in combustion rates relative to the deflagration 
mode. The detonation wave has the added benefit of attaining a higher peak pressure. 

Indeed, the interest in PDE’s arises primarily from advantages that accrue from the significant 
combustion pressure rise that is developed by the detonation process. Conventional rocket engines, for 
example, must obtain all of their compression from the turbopumps, whereas pulse detonation engines 
provide additional compression in the combustor. Thus, PDRE’s are expected to achieve higher specific 
impulse (/ ) than conventional rockets and require smaller turbopumps. The increase in l sp and the decrease 
in turbopump capacity must be traded off against each other to determine the optimal engine choice. Similar 
arguments apply with respect to airbreathing applications. 



A major objective of current development work is to demonstrate that PDE’s can be operated at 
reduced levels of precompression such that flowpath turbomachinery can be greatly reduced. In this case, 
compression, combustion, and thrust-generating processes can be consolidated within a single component, 
yielding a high thrust-to-weight ratio, relative design simplicity, and enhanced reliability. 

The scientific feasibility of the basic PDE concept has been proven in several experimental and 
theoretical studies, as noted in recently published review articles. 1 ’ 2 Significant strides are also being made 
with respect to the development of design strategies and performance optimization. 3 However, there are 
certain developmental issues affecting the application of pulse detonation propulsion systems which have 
yet to be fully resolved. 

Practical detonative combustion engines, for example, require a repetitive cycle of charge induction, 
mixing, initiation/propagation of the detonation wave, and expulsion/scavenging of the combustion product 
gases. Clearly, the performance and power density of such a device depends upon the maximum rate at 
which this cycle can be successfully implemented. The key point here being that accomplishment of these 
tasks requires power that must be supplied by either the engine or delivered from an auxiliary source. 

For instance, the energy needed for direct detonation initiation may be considerable, depending 
upon the specific fuel-oxidizer combination, and the pumping power required for effective scavenging of 
the chamber is a critical design consideration. PDE’s, however, have no shaft power; therefore, an external 
power source or utilization of an alternative energy conversion scheme must be relied on. As such, this 
research investigates the opportunity for direct electric power production using the detonation combustion 
wave to drive a magnetohydrodynamic (MHD) generator integrated into the engine design. 

The author is aware, of course, that a detonation-driven MHD generator can only convert the kinetic 
energy of the combustion products into electrical power. Thus, the enthalpy extraction efficiency of such a 
device will not be comparable with a typical generator configuration utilizing thermal combustion conditions 
and a flow-accelerating nozzle. Nevertheless, the detonation-driven MHD generator may, despite its intrinsic 
efficiency limitations, prove effective for specialized airborne missions. 

A review of the open literature revealed only one paper directly addressing this concept on an 
experimental basis, and it was published by researchers from the Krzhizanovsky Power Institute in Moscow. 4 
Their motive at the time was the development of peak/emergency power stations which would operate on 
natural gas. Because detonation-driven MHD power avoids the need for a high-temperature heat exchanger, 
a large compressor station, and a complex dc-to-ac conversion system, they realized that the intrinsic 
energy conversion inefficiencies associated with this device could be offset by the extremely low capital 
costs of construction. 

These researchers used a natural gas/0 2 -fired detonation tube seeded with a potassium carbonate 
(K 2 C0 3 ) aqueous solution to drive a 16-mm-diameter linear MHD channel subjected to a 0.23 T magnetic 
induction and were able to extract electrical power while running the system at repetition rates approaching 
100 cps for a duration of 2 hr. The channel electrodes were segmented, but the active MHD length was not 
specified. Observed cycle-to-cycle variations in detonation velocity did not exceed 1.5 percent, and the 
wall heat flux accounted for 25 percent of the released combustion energy. The power output peaked at 
—0.1 W e at a load impedance near 10 Q. The maximum effective electrical conductivity measured during 
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these tests was (c) = 3.3 S/m. These test results validated the fundamental technical feasibility of pulse 
detonation-driven MHD power and provided confidence that the concept could be extended for application 
to aerospace propulsion systems. 

The fundamental objective of this line of research was to address developmental concerns which 
affect the prospects for realizing an integrated pulse detonation propulsion and MHD power system. In this 
technical publication (TP), the energy requirements for direct detonation initiation of various fuel-oxidizer 
mixtures as well as the pumping power requirements necessary for effective scavenging of the chamber are 
critically examined. In addition, this TP reports a series of experiments designed to investigate the basic 
engineering performance characteristics of a laboratory-scale pulse-detonation-driven MHD generator. 
This included quantification of ionization properties behind the detonation wave and evaluation of electrical 
power extraction characteristics from a short-length, continuous-electrode Faraday channel. A simple 
electrodynamic model incorporating near-electrode potential losses is developed and is used to correlate 
the experimental work and to extrapolate performance at increased scales. The experimental results reported 
here are only for stoichiometric mixtures of acetylene (C 2 H 2 ) and 0 2 with an atomized spray of cesium 
hydroxide (CsOH) dissolved in alcohol as an ionization seed in the active MHD region. As part of the 
experimental work, an engine stand was developed which allowed direct time-resolved measurements of 
system thrust both with and without MHD interaction. 
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2. OPERATIONAL POWER REQUIREMENTS 


The power required to operate a PDE is dominated by the electrical power consumed in the detonation 
initiation process and by the power needed to transport the propellants. Without loss of generality, assume 
that the pumping power can be delivered using an electric drive. The anticipated power consumption 
requirements for these operational functions are now examined. 

2.1 Direct Detonation Initiation 

Detonation initiation may be accomplished through one of two possible modes: (1) A slow 
deflagration-to-detonation transition (DDT) via an accelerating combustion wave and (2) a fast, direct 
initiation process in which the blast wave from a strong igniter produces the immediate onset of detonation. 
Because the flame travel necessary for DDT transition is on the order of 1 m for the most sensitive fuels 
available, the slow DDT mode does not appear promising for PDE applications. Inquiry is therefore confined 
to the direct initiation mode in which the major parameter of interest is the critical igniter energy. A thorough 
and complete discussion with respect to both modes of detonation initiation may be found in the survey 
paper of Lee.^ 

The most common semiempirical methods for determining the critical detonation energy are based 
on the linear detonation tube method first proposed by Zeldovich et al. in which a spherical detonation is 
initiated by firing a planar detonation wave from a linear tube into a larger volume. 6 Experimentally, it is 
found that the planar wave will transform into a spherical detonation wave if the tube diameter d exceeds a 
certain critical value d c (i.e., d > d c ). If d < d c , the reaction zone will decouple from the shock, and a 
spherical deflagration wave results. 

The critical tube diameter is known as a dynamic parameter of detonation in that it is governed by 
both the chemical reaction rates and the three-dimensional structure of the detonation wave and cannot be 
determined from the equilibrium (static) Chapman-Jouget (C-J) theory. 7 Other dynamic parameters include 
the cell size A, the initiation energy VV ? , the chemical induction time q, the induction length A = at; (where 
u is the particle velocity relative to the shock), and the detonation limits. It turns out that an extensive array 
of experimental evidence indicates a universal linear relationship between d c and A; that is, d c = 13 A. This 
universal correlation is extremely useful in practice when available experimental apparatus permit the 
direct determination of one parameter rather than the other. 

Once A or d c has been experimentally determined, a number of semiempirical point blast initiation 
models have been developed which can be used to estimate the critical initiation energy W s for an unconfined 
spherical detonation wave. Some of these models are described and compared against experimental data in 
the work of Benedick et al. 8 A particularly simple but useful correlation is the “work done” model of Lee 
and Matsui. 9 This model assumes that as the wave emerges from the tube, the energy is delivered to the 
outside volume via the work done by the interface that separates the expanding combustion products from 
the unreacted mixture. A simple expression is thus obtained relating W s to both d c and the C-J states of the 
mixture: 
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w _ *Pcj u cA (1) 

5 24 a cj 

where pQj is the C-J detonation pressure, « CJ is the C-J particle velocity, and a CJ is the speed of sound 
behind the C-J detonation. 

It is also common practice to represent the critical initiation energy in terms of a dimensionless 
number D H formed by the ratio of the critical initiation energy of the fuel-oxidizer mixture of interest to 
that of C 2 H 2 /0 2 . Thus, D H represents a relative measure of detonation sensitivity. 

Another extremely useful parameter is the explosion length R 0 associated with a given source energy: 


* 0 = 



V/O+i) 


ypoj 


( 2 ) 


where; = 0, 1, 2 for planar, cylindrical, and spherical geometries, respectively; Wj is the source energy per 
unit area, per unit length, or just the energy for the three geometries; and p 0 is the initial gas pressure. The 
importance of R 0 is derived from the fact that it is equivalent for all three geometries. 5 Thus, once the 
critical initiation energy is known for the spherical detonation, it is possible to obtain the initiation energy 
for the planar ( W ) and cylindrical (W c ) geometries through the following scaling relationship: 


W, W 

^ w c w p ■ 


(3) 


Table 1 represents the dynamic detonation parameters for a broad range of fuel-oxidizer mixtures 
as deduced through the application of the “work done” model and through the explosion length scaling 
relationship to the experimental data of Matsui and Lee 10 and Beeson et al. 1 1 These data correspond to near 
stoichiometric compositions at 1 atm pressure. 

It is evident from these results that the detonation initiation energy required for H 2 /0 2 mixtures is 
roughly 10 3 times greater than that for C 2 H 2 /0 2 mixtures. This is a minimal requirement and can easily be 
achieved with an auxiliary onboard power supply. For an H 2 /air mixture, however, the detonation initiation 
energy requirement becomes several orders of magnitude higher and the use of auxiliary onboard power 
becomes more problematic. Note that the results for an n-hexane/air mixture may be taken as representative 
of a jet fuel/air mixture and that the electrical power demand for detonation initiation in this case will also 
be significant. In general, the conclusion reached is that the use of practical fuel-oxidizer mixtures in high 
firing rate PDE’s will require a substantial amount of electrical power for direct detonation initiation. This 
power demand can be minimized by using an initiator tube of the smallest feasible diameter (as defined by 
d c ), but it is evident that an integrated power generation system would prove advantageous. 
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Table 1. Dynamic detonation parameters at near stoichiometric composition. 


Mixture 

d c (cm) 

X(mm) 

o H 

W S (J) 

fl 0 (m) 

W p (J/cm 2 ) 

Acetylene/Oxygen 

(C 2 H 2 /0 2 ) 

0.158 

0.122 

1 

3.83 xIO ” 4 

.1.56 xIO " 3 

1.57x10-2 

Ethylene-Oxide/Oxygen 

(C 2 H 4 0/0 2 ) 

0.3 

0.231 

3.1 x 10 ’ 

1 . 2 x 10- 2 

4.91 x10 “ 3 

4.98x10-2 

Ethylene/Oxygen 

(c 2 h 4 /o 2 ) 

0.647 

0.498 

1.9 x10 2 

7.2x10 - 2 

8.92x10- 3 

9.05x10- 2 

Propylene/Oxygen 

(C 3 H 6 /O 2 ) 

0.841 

0.647 

5.3 xIO 2 

2.03x10-’ 

1.26x10- 2 

1.28x10-’ 

Propane/Oxygen 

(C 3 H8/0 2 ) 

1.01 

0.777 

1.5 x 10 3 

5.77x10-’ 

1.78x10-2 

1.82x10"’ 

Ethane/Oxygen 

(WO;,) 

1.1 

0.846 

2.8 xIO 3 

1.07 

2.19x10- 2 

2.23x10-’ 

Methane/Oxygen 

(CH 4 /0 2 ) 

5.28 

4.06 

1.3 xIO 5 

5.07x10-’ 

7.94x10 

8.04x10-’ 

Hydrogen/Oxygen 

<h 2 /o 2 ) 

1.99 

1.53 

4.1 xIO 3 

1.58 

2.5x10-2 

2.53x10-’ 

Acetylene/Air 

(C 2 H 2 /Air) 

8.1 

6.23 

3.4 xIO 5 

1.29 xIO 2 

0.108 

1.11 

Ethylene-Oxide/Air 

(C 2 H 4 0/Air) 

30 

23.1 

2 x 10 7 

7.62 xIO 3 

0.422 

4.28 

Ethylene/Air 

(C 2 H 4 /Air) 

80 

61.5 

3.1 x 10 s 

1.2 x 10 5 

1.06 

10.7 

n-Hexane/Air 

(C 6 H 14 /Air) 

144 

111 

1.8 xIO 9 

7 x10 s 

1.9 

19.4 

Propylene/Air 

(CsHe/Air) 

150 

115 

2 x 10 9 

7.55 x10 s 

1.95 

19.8 

Propane/Air 

(WAir) 

220 

169 

6.6 xIO 8 

2.52 xIO 6 

2.92 

29.6 

Ethane/Air 

(C 2 He/Air) 

280 

215 

1.3 x 10 10 

5.09x10® 

3.69 

37.4 

Methane/Air 

(CH 4 /Air) 

1,020 

785 

5.9 xIO 11 

2.28x10® 

13.1 

133 

Hydrogen/Air 

(H 2 /Air) 

196 

151 

1.1 xIO 10 

4.16x10® 

3.45 

35 


‘All parameters deduced from experimental data of Matsui and Lee 10 and Beeson et al . 11 
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The critical energy of detonation is also dependent upon the rate at which energy is deposited by the 
igniter (i.e., igniter power). In general, the critical energy decreases with a decrease in the duration of 
energy release t ig (i.e., as igniter power increases). This behavior has important implications with respect to 
the optimal design of an igniter circuit for a PDE. 

A simple theoretical model for determining the correlation between critical energy and power required 
for direct initiation has been formulated by Abouseif and Toong. 12 This model considers the flow generated 
by the motion of a constant velocity shock wave in planar, cylindrical, and spherical geometries. The shock 
increases the temperature and pressure of the detonable mixture, and ignition can occur in the shock-heated 
gas after a period corresponding to the induction time. If this blast wave is sufficiently strong, the reaction 
zone is able to couple with the shock wave to form a detonation wave. 

The basis for the model is the observation that a constant velocity shock can be generated by a 
constant velocity piston in the three geometries of interest. 13 Furthermore, Chu demonstrated that such a 
constant velocity piston could be realized in practice by means of an appropriate energy deposition rate. 14 
Abousief and Toong took this source power relationship for generating a constant velocity piston motion 
and hypothesized that the chemical induction time for detonation is equivalent to the duration of piston 
motion for initiating detonation. This hypothesis was verified for the planar case through numerical 
calculations. Following this hypothesis, the necessary power P of the energy source is obtained from Chu’s 
relationship in the form 


P = 


7 

7-1 




where C- - 1, 2k, 4k for j = 0, 1,2 for the three geometries, respectively; p p and u p are the pressure and 
velocity at the piston surface, y is the ratio of specific heats, and r ; is the induction time. The total deposited 
energy is obtained by integrating the power relationship over the induction time and it follows that 


W = — — -i—p uJ +l ti +1 

c y _ lj + l PpP « 


(5) 


The critical energy is time dependent for the cylindrical and spherical geometries. In order to relate 
the critical source energy to a critical source power, it is useful to define an average or effective power as 



7 C j 

7 — 1 j + 1 


t'p P l 


( 6 ) 


Determination of the pressure and velocity of the piston is simple for the planar geometry. In this 
case, the entire flow field bound by the piston and the shock wave is uniform and identical to that just 
downstream of the shock. Thus, p p and u p are equivalent to the conditions just downstream of the shock 
front (i.e., p s and u s as given by the Rankine-Hugoniot relations for a given blast wave strength of Mach 
number M): 15 
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( 7 ) 


p s _ 2 yM 2 +1-7 

po r+i 

T s [27 M 2 + 1 - y|(y - 1 ) M 2 + 2] 
7J) (y + l) 2 M 2 


= M 1 (r-i)A/ 2 +2 

a o (r+i)M 2 


(9) 


Note that specification of the shock Mach number yields a post-shock temperature T s that can be used to 
deduce t t . 

For a one-step, first-order Arrhenius rate expression, the chemical induction time is proportional to 
an exponential function involving the activation energy E a and the temperature T s : 


t i exp 



( 10 ) 


For estimation purposes, a ratio is formed between the induction time for any blast wave strength to the 
induction time t i CJ for a C-J wave such that 


t i= t i,CJ ex PI 


E a( 1 


R 


K T s 


l CJ J 


(ID 


where T CJ is the temperature behind a C-J shock wave. 

This methodology is demonstrated for stochiometric H 2 /0 2 and H 2 /air mixtures since they are of 
great practical interest. The activation energy for H 2 /air is taken as E a = 16.5 kcal/mol, as suggested by 
Shepherd. 16 The activation energy for H 2 /air mixtures was taken as E a = 32 kcal/mol, as deduced by Moen 
et al. 17 The C-J induction times are estimated by fitting available kinetic data. The results of these fits are 
shown in figure 1 where t i CJ ~ 8.8 x 10~ 8 sec for H 2 /0 2 and t t CJ * 2.5 x 10 -5 sec for H 2 /air. 

Predictions for the critical initiation energy as a function of average deposition power are shown in 
figure 2 for a planar geometry. As anticipated, the blast wave strength is seen to increase with power input 
and the critical initiation energy asymptotically approaches a minimum threshold value. It is again clear 
that the electrical power requirements for H 2 /air mixtures can be more than 10 times greater than for H 2 /0 2 
mixtures. 
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Figure 1. Variation of chemical induction time with blast wave strength for stoichiometric 
H 2 /0 2 and H 2 /air mixtures. Based on empirical fit of C-J induction time to available 
kinetic data. 




Figure 2. Predicted critical initiation energy versus average deposition power for planar 
detonations in stoichiometric H 2 /0 2 and H 2 /air mixtures. 


2.2 Chamber Scavenging 

A distinguishing feature of PDE’s is that the blowdown of gas from the chamber during every pulse 
is coupled with a scavenging and recharging process. Such operation is possible only when a propellant- 
pumping mechanism is available. In general, the cylinder pressure will briefly fall below the exhaust pressure 
p e after the blowdown process due to the inertia of the exhaust gases. Once the inlet valve is open, fresh 
mixture will begin to flow into the chamber as soon as the chamber pressure falls below the scavenging 
pressure p r This inflow of fresh mixture continues as long as the inlet valve remains open and the inlet total 
pressure exceeds the pressure in the chamber. 
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It has been suggested that a favorable scavenging effect may be obtained through proper tuning of 
the pressure waves in the exhaust and inlet sections of the engine. That is, the chamber vacuum pulled by 
the blowdown process may be strong enough to eliminate the need for a scavenging compressor such that 
self-aspirating operation becomes feasible. However, it seems unlikely that a design strategy of this type 
can be implemented over a wide range of operating conditions. In all probability, some pumping power 
will be required to achieve effective scavenging and recharging. 

In an ideal scavenging process, the fresh mixture would push the residual chamber gases out of the 
tube without mixing or exchanging heat. This process would continue until the fresh mixture just fills the 
chamber without escaping through the exhaust. In an actual engine, of course, some mixing and heat 
exchange would occur and some of the fresh mixture would be lost through the exhaust. Furthermore, a 
nonreactive buffer gas layer would probably be required between the fresh mixture and residual gases in 
order to prevent preignition. This would represent even further degradation of the scavenging process. 

The effectiveness of the actual scavenging process with respect to the ideal process can be expressed 
in terms of some predefined parameters. For example, the fresh mixture supplied in the ideal scavenging 
process would be that which would just fill the chamber at inlet temperature and exhaust pressure (i.e., 
NV c p s where N is the firing rate, V c is the chamber volume, and p s = p s (T i ,p e ). 

Thus, a scavenging ratio R s may be defined as the ratio of the actual mass supplied to the ideal mass 
supplied: 


R — a 
S ^cPs 


( 12 ) 


A scavenging efficiency e s may also be defined as the ratio of the actual mass of fresh mixture 
retained in the chamber to that for an ideal scavenging process: 


m a 

£S ~ WcPs 


( 13 ) 


T: 


These two parameters may be related to each other through the introduction of a trapping efficiency 


e s = TR s ■ < 14 > 

It is useful at this point to consider an extreme case where the fresh mixture mixes completely with 
the residual gases under the assumption that the residual gases have the same temperature and molecular 
weight as the fresh mixture. For this perfect mixing model, mass conservation may be expressed in the 
following differential form: 


dv - e s dv = V c de s 


( 15 ) 
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where v is the volume of fresh mixture which has flowed into the chamber up to a given instant in time. 
Separating variables and integrating over the scavenging period yields 


ln(l-e 5 ) = 


Z^s 

V, 


(16) 


But R = v JV C by definition and the following working relationship is obtained: 


e=\-e Rs . 


(17) 


Substituting equation (14) yields an additional expression for the variation in trapping efficiency with 
scavenging ratio for the perfect mixing model 


R s 

Figure 3 shows the variation in e s and T with R s for the perfect mixing model. For an ideal scavenging 
process, T = 1 and e s = R s . In practice, e s and R s could be measured and used as a quantitative measure of 
the scavenging effectiveness associated with a particular engine design. 



Figure 3. Variation in scavenging efficiency e s and trapping efficiency T 
with scavenging ratio R s for a perfect mixing model. 
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Although the scavenging efficiency and scavenging ratio are highly important characteristics of a 
PDE, the power required to attain a given scavenging ratio is of equal or greater interest. For the purpose of 
analysis, consider the flow train sketched in figure 4 as an effective representation of an actual engine. 


Scavenging 

Compressor 



Figure 4. Schematic representation of PDE flow train. 

First note that the power to scavenge is a function of the mass flow of fresh charge and the pressure 
at which this charge must be delivered to the inlet port of the chamber. The mass flow is given by 


lh a =A c Ca iPi < P ’ (19) 

where A c is the chamber cross-sectional area (i.e., the characteristic cross-sectional area), C is the engine 
flow coefficient, a t is the speed of sound in the fresh mixture, and (j) is the compressible flow function: 


y + 1 
7 


By substituting equation (19) into equation (12), one obtains the expression 


0 = 


7-1 




\Pij 


A 


V ) 


( 20 ) 




Caj ] 

hi 

Wc) 

yPs J 


0 . 


( 21 ) 


where L c is the chamber length. It is now possible to replace p i /p s with p i lp e since both densities are 
evaluated at the inlet temperature. The result is a relationship for the pressure ratio required for scavenging 
as a function of characteristic engine parameters: 


p i _ R s NL c 1 ( 22 ) 

P e Ca i 0 

Figure 5 shows the variation in scavenging pressure as a function of a characteristic engine velocity 
assuming air as a working fluid. Note that the pressure differential required for effective scavenging can 
become significantly large under some circumstances, implying the need for high pumping power. 
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Figure 5. Variation in scavenging pressure ratio as a function of characteristic 
engine velocity. Assumes air as the working fluid. 


For a given required scavenging pressure and inlet pressure pq, it is possible to estimate the 
necessary pumping power. For adiabatic, isentropic compression of an ideal gas, the shaft work w ca may be 
expressed as 

»ca=W<?c ■ <23) 

where C p 0 is the constant pressure specific heat and Y c is a pressure ratio function for the compressor: 

r - 1 




V 

yPoy 


-1 . 


(24) 


Therefore, the power required for adiabatic, isentropic compression of an idea gas is given by 


P ca = m i W ca = lh i C P 0 T 0 Y c 


(25) 


The actual shaft work w may be expressed as 


w c~ C po{ T i~ T o) ’ 

and the actual power required for compression is given by 


(26) 


P c = rh i w c = m i C p0 {T i -T 0 ) 


(27) 
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At this point it is beneficial to define a compression efficiency as the ratio of the adiabatic shaft 
power to the actual shaft power: 

p T V 

T] = — = -0-C— , (28) 

C Pc Ti-To 

such that 


and 



ril i C pQ T 0 Y c 

Pc 


(29) 


T; 


1+-^ 

P. 


C J 


(30) 


An energy density parameter may now be formed by substituting equation (12) into equation (29): 


jc _ (3D 

N Ps V C R S Pc 

Figure 6 shows the pumping power per unit mass flow as a function of the scavenging pressure ratio p^Pq, 
assuming air as the working fluid with r\ c = 0.9. Note that equation (3 1) is a useful scaling relation in terms 
of scavenging parameters and demonstrates the level of pumping power that might be required in a practical 
engine. 



Figure 6. Variation in pumping power per unit mass flow as a 
function of the scavenging pressure ratio. Assumes 
air as working fluid with rj c = 0.9. 
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3. THEORETICAL AND EXPERIMENTAL PERFORMANCE 


3.1 Experimental Apparatus 

The general experimental arrangement is shown schematically in figure 7. The laboratory-scale 
detonation tube is closed at one end and various plasma diagnostic and MHD channels can be attached to 
the open end using an area transition adapter. A permanent magnet assembly and an electromagnet were 
custom fabricated for establishing MHD interaction in the channel. All of the experiments were conducted 
in a single-shot mode where the system is evacuated, charged to atmospheric pressure with premixed 
combustibles, and ultimately detonated using a high-voltage spark discharge near the closed end of the 
tube. During this process, a thin membrane separates the chamber volume from the surrounding atmosphere 
until it is ruptured by the detonation wave. Ionization seed is introduced using an atomizing injector that is 
located near the open end of the detonation tube and oriented in such a way as to direct the seed spray along 
the length of the attached channel. The entire tube/channel assembly is suspended on rollers and preloaded 
against a piezoelectric-type dynamic load cell. Dynamic pressure transducers were installed in the tube to 
measure detonation wave strength and to provide time-of-flight measurement of detonation wave speed. 
The design/configuration of the individual components making up the experimental apparatus are discussed 
in detail in the following subsections. 



Diaphragm 


Figure 7. General schematic of laboratory-scale pulse detonation MHD apparatus. 
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3.1.1 Detonation TYibe 


The detonation tube is constructed from a 1-m-long section of schedule-40 steel pipe having a 
nominal inside diameter of 2.54 cm. The detailed configuration of this tube is shown in figure 8. The head 
of the detonation tube is sealed with a pipe cap in which fitting connectors have been installed for evacuating 
the tube and introducing premixed combustibles. The open end is flanged for attachment of an area transition 
adapter and diagnostic/MHD channels. 



Figure 8. Schematic of the 1-m-long, laboratory-scale pulse detonation tube with 2.54-cm bore. 


A pair of high-voltage igniter electrodes were fabricated from copper (Cu) pins with tungsten (W) 
tip inserts and located near the closed end of the tube. The tip of the igniter cathode was ground to a point 
to increase the local electric field intensity and thereby reduce the required breakdown voltage across the 
3-mm gap. The discharge energy for directly initiating the combustion detonation wave is created by 
discharging the energy stored in an inductive-resistive-capacitive (LRC) circuit that is closed using a solenoid- 
operated mechanical switch. The oscillating electric field that drives the gas discharge was observed to 
decay completely within 100 //sec. The circuit discharges a 2-//F energy storage capacitor which has been 
precharged to 10 kV, yielding an ignition energy source of magnitude W ign = */2 CV 2 ~ 100 J. This value 
exceeds the minimum ignition energy necessary for direct initiation of C 2 H 2 /0 2 mixtures by several orders 
of magnitude. 10 This level of ignition energy provides assurance that true detonation conditions are attained. 

Two dynamic pressure transducers are mounted in the tube at roughly equal distances from either 
end such that the separation distance between them is exactly 0.25 m. These sensors provide an indication 
of detonation wave strength and are used for time-of-flight measurements of detonation wave speed. For a 
data sampling rate of 200 kHz, an uncertainty of ±5 //sec is encountered in the time-of-flight measurement. 
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3.1.2 Atomizing Seed Injector 

An electrically activated atomizing fuel injector from an automotive engine was adapted for injection 
of ionization seed. The injector was mounted to the detonation tube and oriented so that the spray cone was 
directed downstream into the channel region. A 30-percent solution (by mass) of CsOH in a methanol 
carrier was used for seeding. In conducting each test, the pressurized injector for 50 msec was first activated 
and then the ignition trigger was delayed by 100 msec. Based on a calibrated volumetric injection rate of 3 
cm 3 / sec, it was estimated that 0.15 cm 3 of seed solution was injected per shot. It was impossible to ascertain 
the actual seed fraction in the channel, however, due to poorly defined atomization characteristics, attachment 
of the seed material to the wall, and uncertainty in the seed distribution. 

3.1.3 Plasma Diagnostic Channel 

A plasma diagnostic channel was also designed and fabricated to quantify the electrical properties 
of the ionized gas behind the detonation front. This short length, axisymmetric diagnostic channel, which 
has the same internal diameter as the detonation tube, is shown in figure 9. The main body is machined 
from G-l 1 phenolic and is fitted with a set of 1.6-mm copper electrode segments to which a preset electric 
potential can be applied. These electrodes are mounted flush with the wall and maintain a contact angle of 
60° with the flow. The channel is also equipped with two tungsten wire electrostatic probes for measuring 
the local plasma potential in the core. They are located 5 mm above the anode and cathode surfaces, 
respectively. Together, the discharge electrodes and electrostatic probes form a four-electrode electrical 
conductivity probe. 18 A tap for a dynamic pressure transducer is located at the same axial location. In 
addition, a sapphire window port was integrated into the design so that the gas discharge process and 
detonation wave passage could be observed optically. The diagnostic channel was fabricated for direct 
mating with the detonation tube. 



Figure 9. Schematic of windowed axisymmetric plasma diagnostic channel. Note that the wall- 
mounted discharge electrodes can be combined with the electrostatic probes in the 
core to form a four-electrode electrical conductivity probe. 
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3.1.4 Laboratory Magnets 


A permanent magnet assembly (0.6 T) and a laboratory electromagnet (0.95 T) were custom 
fabricated to support testing efforts. The permanent magnet is composed of individual segments as shown 
in figure 10. Each segment is 2 in. wide and contains eight 2 x 2 x 1 in. low-grade, rare Earth magnetic 
blocks ( BH max = 35 MGOe) and four 2 x 2 x */2 in. high-grade, rare Earth magnetic blocks (BH max = 45 
MGOe). A total of six segments combine to form a 30-cm-long, permanent magnet assembly that provides 
a uniform magnetic induction of 0.6 T over a 3.8-cm air gap. The laboratory electromagnet is shown 
schematically in figure 11. The bore for this magnet is 46 cm long and has a 5-cm air gap. The maximum 
magnetic induction in the bore is 0.95 T at full design current. 




in 


1010 Steel Pole Pieces 
Permanent Magnet (35 MGOe) 
Permanent Magnet (45 MGOe) 
Spacers 


Figure 10. Cross section of permanent magnet segment. Six segments are combined to 
form a 30-cm-long assembly having a uniform magnetic induction of 0.6 T 
over a 3.8-cm air gap. 


A- — i 




Figure 11. Schematic of the laboratory electromagnet. The magnet bore is 46 cm long 
with an air gap of 5 cm. The magnetic induction in the bore is 0.95 T at full 
design current. 
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3.1.5 Faraday MHD Channel 


A continuous-electrode Faraday MHD channel was also designed and fabricated to fit the available 
laboratory magnets. This constant-area channel, shown in figure 12, has an active MHD length of 30 cm. 
The internal flow passage is rectangular in shape with a 2.54-cm height between electrodes and a 2-cm 
width between sidewalls. The main body of the channel is fabricated from G-l 1 phenolic, as is the area 
transition adapter connecting the tube to the channel. The copper alloy electrodes include pressure taps at 
four locations along the active channel length. 



Figure 12. Schematic of the 0.3-m-long, continuous-electrode Faraday MHD 
channel. Active region is 2.54 cm in height and 2 cm in width. 


3.1.6 Experiment Control and Signal Acquisition 

A PC-class computer with appropriate I/O hardware was used to sequence and control the experiments 
and acquire four channels of high-speed data at a sampling rate of 200 kHz per channel. A 50-MHz, two- 
channel digital oscilloscope was also utilized to capture fast signal transients. High-speed diagnostic 
measurements available for physical interpretation include the following: Time-resolved thrust profiles as 
a means of deducing I performance; static pressure fluctuations at the wall to define the speed and strength 
of the combustion detonation wave; radiative emission through the window port of the diagnostic channel 
for additional definition of detonation wave morphology and structure; electrostatic probe voltages for 
measuring local plasma potential relative to the cathode; anode-to-cathode voltages and currents in the 
plasma diagnostic channel circuit as a means of deducing the electrical properties of the ionized gas; and 
load voltages and currents in the MHD channel circuit as a means of deducing power extraction 
characteristics. 
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3.2 Electrodynamic Model 


A simple, first-order electrodynamic model was developed as an aid to the interpretation of 
experimental results and as a basis for developing generalized scaling laws. In this model, it is assumed that 
a uniform slug of burned ionized gas moves with velocity u through a channel of height h that is subjected 
to an applied magnetic induction B. Wall electrodes directly couple MHD interaction in the flowing plasma 
with an LRC active loading circuit. In general, the slug thickness S may be less than the active channel 
length L. A uniformly distributed current density j in the plasma slug is assumed, and Hall effects are 
neglected due to the very high-pressure levels experienced behind the combustion detonation wave. 
Furthermore, the magnetic Reynolds number is assumed to be much less than 1, and flux compression 
effects were neglected. The electrodynamic model and potential level diagram are shown in figure 13. The 
active loading circuit includes the possibility of precharging the capacitor to an applied voltage v 0 as a 
means of simulating a higher effective magnetic induction, B e ^ 



E/nd=e//J = uxB 


Eapp = Wh 


E = (//?+ //?/)//) 


E d= iR<j/h 


Figure 13. LRC active loading circuit and electric potential level diagram 
for idealized electrodynamic model. 


The effective induced electric field E e ^ in the channel is the sum of the induced Faraday field 
associated with the applied magnetic induction, E ind = | u x B | = dh, and the applied field, E = 
Thus, the effective induced potential in the channel may be written as 

£ e ff =£ + v 0 = uBh + v Q = uB e ffh , ( 32 ) 
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where 


• 

Requiring that the current have a positive magnitude in the direction of u x B, Kirchoff’s voltage 
law implies that the effective induced electromagnetic force (emf) £ e ff= uB e ^h be balanced by voltage 
drops resulting from the external load R, internal resistance of the plasma /?•, and a lump potential loss v d = 
iR d : 


£ e jj = uBgffh. = iR + iR[ + i R d = V + jh / <7 + v d , (34) 

where the current density,./ = HA, and the electrical conductivity, cr= h/(AR ,), have been introduced. The 
lump loss term v d is introduced to account for plasma nonuniformities, near-electrode potential drops, 
emission losses, current concentrations, as well as possible current leakage and shortings. 19 It is assumed 
that v 0 « constant during passage of the plasma slug, and because the inductance in the loading circuit is 
extremely low, inductive energy storage terms have been neglected. Thus, equation (34) may be further 
manipulated to obtain 


iiBgjj — V/h + j/<7 + v d /h — E + j/(T + E d , (35) 

where i is the load current and V is the terminal voltage measured across the load. Solving equation (35) for 
j yields the generalized Ohm’s law for the generator: 

j = <\«B eJf -E-E d \ . ( 36 ) 

Following a modeling strategy analogous to Wu, 19 a dimensionless voltage drop A is defined as 

A = v d juB e ffh . ( 3? ) 

This engineering parameter represents a direct measurement of the voltage drop losses relative to the 
effective emf of the generator. 

For an MHD generator, the load factor K is defined as 


electric power delivered to load i 2 R 
useful electric power generated i 2 R + i 2 R i 

and elimination of v d for A yields the relation 


iR V _ E 

iR + iR i uBgjjh — v d uB e ^{\- A) 


(38) 


( 39 ) 


21 



By substituting equations (37) and (39) into the generalized Ohm’s law, an expression is obtained for the 
current density in the generator: 

j = CuB eff {\-A)(l-K) . ( 4 °) 

In addition, an expression for the generator power density follows from the basic definition 
tP = j • E and from use of equations (39) and (40) to eliminate E and j: 

= j . E = Ou 2 B 2 eff {\ - A) 2 £(l - K) . ( 41 > 

The maximum power density occurs when dtP/dK = 0 such that K = ! /2 for optimum performance. It should 
be clear from equation (41) that when the induced emf is limited due to small channel height and/or low 
magnetic induction, the lump potential loss v d can have a profound affect on generator performance. 

The Lorentz body force acting on the plasma slug may be defined in terms of the charged particle 
drift velocity w: 


f = nq(y/ x B) = j x B ( 4 ^) 

and substitution of equation (40) into equation (42) yields the actual Lorentz body force experienced by the 
plasma slug: 


/ = | j x B| = GuB eff B{ 1 - zl)(l - K) < j x B 


'eff 


(43) 


Thus, for the active loading circuit,/ B e ^B assuming all other factors remain constant. However, for true 
similarity, the Lorentz body force should vary in proportion to B 2 e not BB e jj, and the MHD interaction 
level is too small by the factor B e ^/B = l+v^(uBh). 

More directly stated, the magnetic interaction parameter S (work done by Lorentz force/inertial 
energy of flow field) based on the use of an effective magnetic induction is less than that required for 
proper simulation of the push work; that is. 


S ~ R m R h ~ 


_ (J 0 L 0 B eff B 0 


P 0 W 0 


P 0 M 0 


= 5. 


required 


(44) 


where R m (Lorentz force/magnetic body force) is the magnetic Reynolds number, R h (magnetic body force/ 
inertial force) is the magnetic force number, and the subscript zero refers to reference conditions. The 
actual rate at which the ionized gas does push work in the actively loaded channel is 
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? p .fu = m 2 B eff B(\-A)(l-K) . 

When C=0 and there is no applied voltage, B e ff= B and the passive load model is recovered. 



3.3 Experimental Results 

3.3.1 Detonation Characteristics 

As a prelude to operating the detonation tube, detonation calculations were carried out using the 
NASA SP-273 chemical equilibrium code. 20 The stoichiometric mixture ratio was based on the combustion 
products being either carbon dioxide or water, and the initial static pressure and temperature in the unbumed 
gas waspj=l atm and Tj=300 K. Assuming C-J detonation wave structure, this program computes, among 
other parameters, the pressure/temperature ratio across the detonation front, the detonation wave velocity 
D relative to the duct, and the speed of sound of the burned gas a 2 relative to the detonation wave. For C-J 
conditions, the burned gas velocity is sonic relative to the detonation wave (i.e., M 2 =1), and 
11 2 = a i = ^(Y2 r 2 t 2 )- Then > ^ is possible to deduce the burned gas velocity relative to the tube using the 
relation u b = D—u 2 = D—a 2 where u b corresponds to the slug velocity in our electrodynamic model. The 
major results, summarized in table 2, provide a baseline performance benchmark for our experiments. 


Table 2. Summary of calculated denotation wave properties. 


D (m/sec) 

2,425 

32 (m/sec) 

1,317 

Ub (m/sec) 

1,100 

P 2 /P 1 

33.64 

m 

14.04 


The computed detonation wave speed is in good agreement with widely accepted experimental 
measurements. 21 A static pressure p 2 = 33.64 atm (495 psia) and a static temperature T 2 = 4,200 K behind 
the detonation wave based on initial conditions in the tube are expected. 

Experiments were then conducted in our apparatus to confirm the attainment of true detonation 
conditions. As an example, representative pressure waveforms acquired from the dynamic pressure 
transducers during a typical shot are presented in figure 14. The slightly elevated pressure levels near the 
ignition source are attributed to overdriven detonation as caused by use of an extremely high ignition 
energy (=100 J). Once the detonation wave has traversed two-thirds the length of the tube, however, C-J 
conditions are fully attained, and the pressure rise across the detonation front is found to be in close agreement 
with the baseline theoretical prediction. As further evidence, note that the cross correlation of these waveforms 
yields a time delay of 105 psec and a detonation wave speed relative to the tube of 2,380 m/sec. This 
quantity is also in good agreement with the baseline calculations. Repeated experiments showed little 
variation in the preceding results, providing confidence that true detonation conditions were being reliably 
achieved. 
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Figure 14. Representative measured pressure waveforms in the detonation 
tube during a typical shot. 


3.3.2 Plasma Diagnostics 

With the plasma diagnostic channel attached to the end of the detonation tube, further experiments 
were conducted to quantify the ionization properties of the burned gas following the detonation front. Of 
utmost interest were (1) the level of electrical conductivity achieved, (2) the magnitude of near-electrode 
potential losses, and (3) the definition of structural features. 

In these experiments, an electric breakdown potential of 100 V was applied between the wall 
electrodes during passage of the detonation wave. The resulting discharge of current through the ionized 
gas allowed deduction of bulk electrical properties from the measured electrostatic probe voltage waveforms 
and the anode-to-cathode current waveforms. Although highly intrusive, this technique provides an effective 
measurement of major electrical parameters of engineering interest. 

The single, most important parameter indicating ionization quality is the electrical conductivity. In 
this work, two intrusive probing methods were utilized for determining an effective bulk value based on 
Ohm’s law, j = oE. They are classified as two- and four-electrode probes, respectively. 18 Application of 
these probing methods can lead to dubious results, however, due to uncertainty in the current distribution 
within the ionized gas. This uncertainty in current distribution makes the accurate assessment of current 
density and derived quantities such as electrical conductivity problematic. 

Reliable deduction of these electrical parameters requires careful consideration of electrode geometry 
and probe configuration. When using point electrodes, for example, fundamental theoretical considerations 
and a vast array of empirical evidence has demonstrated that the current tends to uniformly diffuse into a 
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spherical distribution having a diameter equal to the anode-to-cathode gap in the ionized gas. 22 Based on 
symmetry arguments for the thin electrode segments used in the diagnostic channel, it is assumed, therefore, 
that the current will uniformly diffuse into a circular cylinder having a diameter equal to the channel height 
ft and a length equal to the electrode width w. This geometry is illustrated in figure 15. 


Electrode 



t 


Figure 15. Illustration of the assumed current distribution between discharge 
electrodes in the plasma diagnostic channel. 


The electric field between the electrodes can be modeled using two-dimensional line sources where 
the electric potential 0 at any point (x,y) from the centerline of the electrodes and the centerline of the 
channel is obtained by summing the contribution due to each line source: 


0 = Aln^jc 2 + (h/2 - y) 2 - A In -^/jc 2 + (h/2 + y) 2 . ( 46 ) 

The constant A represents the source strength and can be evaluated by imposing the known voltages on 
the top and bottom electrodes. The electric field vector E can be determined from the relation E = -grad 0. 
The equipotential surfaces within the assumed cylindrical current volume are illustrated schematically in 
figure 15. 


3.3.2.1 Two-Electrode Probe. In the two-electrode method, consider first a cross-section of the 
assumed cylindrical current volume and define the current density j through an equipotential surface A es as 


j = 



(47) 
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where i is the measured total current and B, has been introduced as an efficiency coefficient of order unity. 
As a simplification to the analysis, approximate the area of the equipotential surfaces as the product of the 
cylinder chord c and the electrode width vv: 


A es ~ wc = wl^fh^-p , 


(48) 


where £ is the radial distance from the electrode surface to the chord, and A es = wh when £ = hi 2. Upon 
substituting Ohm’s law, an expression is obtained for the local electric field in the form 


i 1 

w£<j 2 ^jhC~C 2 


(49) 


Then, integrate for the applied voltage across the channel, using equation (49): 


A/2 

v app = 2 j Ed C 


0 


■ h I 2 Ar 

j f 

0 


(50) 


Evaluation of this integral yields the final desired expression for electrical conductivity in terms of the 
primary measurement parameters: 


K i 

2 v wB 
app ’ 


(51) 


3.3.2.2 Four-Electrode Probe. Unfortunately, the two-electrode probe tends to suffer from near- 
electrode potential drops and nonohmic contact resistances, and a four-electrode probe is usually 
recommended. In the four-electrode probe, a pair of electrostatic probes are introduced between the anode 
and cathode to measure the voltage drop v di j over a distance d in the core of the flow. With this configuration,, 
the parasitic resistance problems are minimized. 

The critical assumption is the definition of the cross-sectional area for computing current density in 
the core. Here, the cross-section of the current volume cylinder at the centerline is taken as the effective 
area such that A e ^ = wh. Thus, the effective current density in the core Uo- takes the form 


i i 


(52) 


The expression for electrical conductivity using the four-electrode configuration follows from a direct 
application of Ohm’s law: 


J jff _ [ d 

E wh % v dif 


(53) 
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3.3.2.3 Ionization Quality/Structure. A seven- waveform ensemble average of the effective current 
density in the diagnostic channel, as computed from equation (52), is shown in figure 16. In this work, £ is 
taken to be unity. These measurements demonstrate the emergence of primary and secondary wave structures. 
The primary wave is associated with the passage of the burned gas behind the detonation front, but the 
physical mechanism responsible for the trailing secondary wave is not clear. The emergence of a conductive 
secondary wave structure has been observed in previous detonation experiments with oxy-acetylene 
mixtures. 4 ’ 23-25 Generally speaking, these experiments indicate a secondary wave structure which may be 
associated with delayed chemical reactions behind the detonation front. As such, the induction period 
never exceeded more than a few hundred microseconds. In figure 16, however, the time delay between the 
primary and secondary waves is rather long (=1.5 msec), and the appearance of the secondary waveform 
shows a temporal correlation with the extreme tail of the measured pressure waveform. Therefore, the 
secondary waveform observed here is more probably associated with flow reversal in the channel. The 
peak current density in the primary and secondary waves was roughly 2 and 0.5 A/cm 2 , respectively. 



Figure 16. Seven-waveform ensemble average of the effective current density 
in the plasma diagnostic channel. Note the delayed appearance 
of a secondary wave structure. 


Normalized seven-waveform ensemble averages of static pressure and radiative emission in the 
diagnostic channel are shown in figure 17. In this case, the hump riding on the trailing edge of the pressure 
waveform is indicative of a secondary compression wave similar to that observed by Cher and Kistiakowsky. 24 
Radiative emission integrated over the optical region is very constant during passage of the primary 
conductive wave, seemingly suggestive of fairly uniform thermal conditions; however, the measured electrical 
conductivity, presented below, does not corroborate this interpretation. 

A seven-waveform ensemble average of the differential electrostatic probe voltage v^yis shown in 
figure 18. Only data corresponding to passage of the primary wave were captured. Note that the voltage 
drop between probes was =40 V during this event. By assuming v d = v app - v di p it was possible to estimate 
the dimensionless voltage drop in the channel using equation (37) where B y = v app /uh. This estimation 
yielded A = 0.6. 
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Figure 17. Seven- waveform ensemble average of static pressure and radiative 
emission in the plasma diagnostic channel. Waveforms have been 
normalized using peak values in the detonation front. Note the 
emergence of a receding Cher-Kistiakowsky secondary compression 
wave. 



Time (msec) 


Figure 18. Seven-waveform ensemble average of the electrostatic probe 
voltage differential in the plasma diagnostic channel. 


Ensemble-averaged, electrical conductivity waveforms are shown in figure 19 based on two- and 
four-electrode probing techniques. Agreement between the two methods is very good. The peak electrical 
conductivity in the primary wave is seen to be ~6 S/m. For comparative purposes, calculations were carried 
out for the electrical conductivity and electron number density in the burned gas as a function of the seed 
mole fraction. This was accomplished using a modified version of NASA SP-273 chemical equilibrium 
code in which the plasma electrical transport properties are computed according to Frost. 26 The results are 
shown in figure 20. 
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Figure 19. Seven- waveform ensemble average of the electrical conductivity in the plasma 
diagnostic channel based on both two- and four-electrode probing techniques. 
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Figure 20. Computed electrical conductivity and electron number density behind the 
detonation front based on a modified version of the NASA SP-273 chemical 
equilibrium code where plasma electrical transport properties are computed 
according to Frost. 

3.3.3 MHD Power Extraction 



MHD power extraction experiments were conducted with the short-length, continuous-electrode 
Faraday channel attached to the detonation tube and placed in the bore of the magnet assemblies. In these 
experiments, the channel was connected to the active loading circuit shown in figure 1 3 in order to simulate 
higher effective magnetic induction. Our principal objectives were as follows: (1) To demonstrate pulse 
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detonation-driven MHD energy conversion, (2) to map the load line characteristics of the generator, and (3) 
to estimate energy density characteristics for engineering scaling purposes. 

3.3.3.1 Open-Circuit Voltage/Burned Gas Velocity. The first task was to acquire open-circuit 
characteristics of the generator by applying a very high load impedance with the capacitor removed from 
the circuit. Figure 21 shows a typical terminal voltage waveform under these conditions, where the waveform 
represents a direct measurement of the peak Faraday potential in the generator, V oc = uBh ~ 10 V. Obviously, 
the small scale of the device greatly inhibits performance in that the near-electrode potential drop is of the 
same order as the Faraday potential. Dividing this waveform by the known magnetic induction (0.6 T) and 
the channel height (0.0254 m), one can deduce an effective burned gas velocity relative to the tube of 
<u b > = 660 m/sec. This value can be compared with the theoretical equilibrium prediction specified in 
table 1 (1,1 00 m/sec). The effective burned gas velocity is lower because the velocity behind the detonation 
wave must decay to zero at the closed end of the tube. 



Time (msec) 

Figure 21. Generator terminal voltage waveform under open-circuit loading 
conditions. Division of this waveform by Bh yields an effective 
burned gas velocity waveform in the channel. 


3.3.3.2 Load-Line Characteristics. The load-line characteristics of the generator were mapped 
out for six values of effective magnetic induction v 0 : 0.6, 0.95, 1.5, 2.4, 3.3, and 4.2 T. For each effective 
magnetic induction, current waveforms for the following set of load impedance values were measured: 0, 
0.0175, 0.098, 1, 5, 10, and 100 £2. Representative results are shown in figure 22 for B eff = 1.5 T. As 
expected, the peak current attains its maximum value under short-circuit conditions and falls with increasing 
load impedance. Also note the emergence of a delayed, relatively weak secondary pulse which is associated 
with the secondary gas-dynamic wave identified in figure 17. 
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Figure 22. Measured current waveforms for various load impedances with B e ^ = 1.5 T. 

Peak current scales inversely with load impedance. 

Generator load lines were constructed based on measured peak values for current and load voltage. 
The results, normalized with respect to the short-circuit current and the open-circuit voltage, are shown in 
figure 23. The characteristic shape exhibits linear features away from the extreme loading conditions. The 
cause of the highly nonlinear behavior near short-circuit conditions is still not understood. 



Figure 23. Measured load-line characteristics of the pulse detonation MHD generator. 
Solid symbols denote true simulation conditions (i.e., B = B e jj). 


31 



3.3.3.3 Energy Density Scaling Parameter. From an engineering point of view, it is advantageous 
and useful to express the power extraction characteristics of the generator in terms of a performance- 
related scaling parameter. Once quantitative estimates are developed for this parameter, it is possible to 
project device performance at practical scales. Although it is difficult to extrapolate the performance of 
MHD devices from laboratory scales, where boundary layer losses dominate, to a scale suitable for practical 
application, where core flow processes dominate, this approach is adopted as a conservative estimate of 
scaling effects. Under virtually all circumstances, this approach will tend to underpredict performance and 
will usually represent a worst-case scenario. The prediction can therefore be treated as an effective lower 
bound on anticipated performance. 

Since electrical power output is the performance parameter of primary interest in MHD generators, 
development is begun with the following definition: 


QL.u21L.aM 

dt dx 


(54) 


where W represents electrical work in Joules and the power density fP was obtained in equation (41). The 
electrical work per unit channel length per detonation pulse can therefore be written as 


dW _ Pwh5 = Wp (55) 

dx u L 

Here, L is the active length of the MHD channel and W D represents the electrical energy transferred per 
detonation pulse. Further note that W D can be empirically evaluated by integrating over the measured 
output waveform of the generator: 


W D = J i 2 R dt + Jljijdt~ J i 2 R dt , (56) 

where the inductance „£is very small in our circuit and can be neglected. 

Equation (55) can now be manipulated to obtain the desired energy density scaling parameter in the 

form 


Wp _ 28 

whL u 


(57) 


This parameter has units of J/m 3 and represents the electrical energy produced per unit channel length per 
unit channel height per unit channel width per detonation. By eliminating fP using equation (41), one 
arrives at the final working relation: 


*D_ 

whL 


(1-A) -2 = SuB 2 K(l-K) 


(58) 
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Note that the parameters on the right-hand side of equation (58) are assumed to be independent of 
device scale. Thus, once the energy density W D /whL is empirically evaluated and one is able to make 
appropriate assumptions concerning the scaling of A, it becomes possible to extrapolate performance at 
increased scales. 

A comparison of measured and theoretical predictions for the energy density scaling parameter is 
presented in figure 24 as a function of load impedance. The empirical value was obtained by integrating 
equation (56) using measured current waveforms and dividing by the known active channel volume. The 
theoretical predictions were computed from equation (58) using the following device parameters: A = 0.6, 
5 = 0.3, <7=6 S/m, u = 1,100 m/sec, B = B e jj= 0.6 + v 0 /uh T, and h = 0.0254 m. Agreement is reasonably 
good given the uncertainty in parameter values. The measured peak electrical energy density ranged from 
10 to 10 3 J/m 3 when the effective magnetic induction was varied from 0.6 to 4.2 T. Clearly, the performance 
of this MHD device is limited by scale and near-electrode potential losses. 



Load Impedence (a) 

Figure 24. Comparison of measured and theoretical predictions for the energy density scaling 
parameter as a function of load impedance. The theoretical predictions are shown 
as curves. Device parameters used in the calculations were as follows: A = 0.6, <5= 

0.3, <7=6 S/m, u = 1,100 m/sec, B e jj = 0.6 + Vq/uIi T, h = 0.0254 m. Optimum 
power extraction occurs when K = 1/2. Solid symbols denote true simulation 
conditions (i.e., B = B e ^). 

The results also show peak power extraction shifting from a load impedance of 10 Q, when 
v 0 = 0 V, to 5 Q, when v 0 = 100 V. This shift is associated with decreasing internal plasma resistance as the 
effective magnetic induction and Joule dissipation increase. That is, optimum electrical power extraction 
always occurs when the load impedance and the internal impedance of the generator are matched. 

At this point, recall that the magnetic interaction parameter S based on an effective magnetic induction 
is less than that required for accurate simulation of the MHD interaction level. That is, the push work done 
by the working medium in these experiments is much lower than would be required if the magnetic induction 
was equivalent to the applied field. In order to justify the extrapolation of these results to the practical case 
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where B = B e g, a forced assumption is made that the enthalpy extraction achievable in the nonsimulated 
case will be high enough to provide the required push work. It is also necessary to assume weak coupling 
between MHD interaction level and burned gas velocity. 

3.3.3.4 Extrapolation of Performance to Practical Scales. The performance of small MHD 
devices tends to suffer from severe boundary layer losses, as evidenced by our experimental results, and it 
is advantageous to use the largest feasible scale for the desired application. Along these lines, the energy 
density scaling relationship defined by equation (58) makes it possible to extrapolate device performance 
to practical scales from the baseline laboratory scale experimental results. As an example, consider a pulse 
detonation MHD generator that is similar in all respects to the laboratory device with the exception of 
channel size and applied magnetic field. Choosing values consistent with practical implementation of this 
concept, a channel height and width of 20 cm, an active interaction length of 1 m, and an applied magnetic 
induction of 2.4 T was selected. 

As a preliminary task to quantifying performance scaling with size, it is first necessary to estimate 
the magnitude of near-electrode potential losses in a practical device. Past experience with MHD devices, 
for instance, indicates that the near-electrode voltage drop can be reduced considerably below that 
encountered in the present device, which was operated only in single-shot mode and therefore suffered 
from severe cold wall effects. In principle, A can be reduced immensely by allowing the electrodes to 
operate at high temperature extremes. Based on past experience, being able to limit the voltage drop to no 
more than 10 percent of the induced field is anticipated. Holding the right-hand side of equation (58) 
constant, an energy density gain above the baseline experimental results in the following estimate: 


("W-M-L oj (l-O-l) 2 ,, (59) 

o.6 d - 0 - 6 ) 2 

For v 0 = 50 V (i.e., B e ^ = 2.4 T), figure 24 yields (W D /whL) A _ 0 6 - 300 J/m 3 at optimal loading conditions. 
Correcting for the anticipated reduction in A using equation (59) results in ( W D /whL) A=0A = (5.1)(300) = 
1,500 J/m 3 per detonation. 

The performance gain associated with increased size follows from the multiplication of the corrected 
energy density with the active channel volume. That is, (W D ) A = 0 A = ( whL)(W D /\vhL) A _ 0>1 = 
(0.2)(0.2)(1)(1,500) = 60 J per detonation. Assuming a practical repetition rate /= 100 cycles per second, a 
continuous power output of magnitude P = f(W D ) A _ 01 = (100)(60) = 6,000 W e is deduced. 

Owing to the poor seeding characteristics of the present device, the preceding value represents a 
very conservative estimate of electrical power generation potential. For instance, the theoretical calculations 
shown in figure 20 indicate that electrical conductivity behind the detonation front can be increased as 
much as 5 times above the level measured in the device. This would imply a direct fivefold increase in 
electrical power generation performance. Optimization of seeding characteristics will require careful tailoring 
of injector configurations and atomization processes so that a uniform seed distribution is obtained along 
the channel length. 
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3.3.4 Thrust Characteristics 


The dynamic load cell was also installed on the test stand and a series of shots were performed, with 
and without MHD interaction, to establish baseline thrust performance. Due to the low conductivity level 
experienced in these experiments, MHD interaction and enthalpy extraction were both very low, and the 
measured thrust profdes were found to be independent of MHD effects. 

An ensemble-average profile of time-resolved engine thrust measurements is shown in figure 25. In 
general, the thrust-producing process may be broken into two phases: (I) A very short period of approximately 
constant thrust during which the detonation wave propagates from the closed end of the tube through the 
stationary mixture to the end of the chamber, and (II) a long period of exponentially decaying thrust associated 
with chamber blowdown. 




Figure 25. Ensemble-averaged thrust profile without MHD interaction and fitted model profile. 


The timescale for phase I is readily determined from the ratio of the chamber length L c to the 
detonation wave speed D (Tj « LJD). The thrust produced during this period is due solely to the pressure 
acting on the closed end of the tube since there is no chamber exhaust until the detonation wave exits the 
chamber. Due to the Taylor wave, 27 a centered expansion occurs between the detonation wave and the 
closed end of the tube, where the velocity must go to zero. Therefore, the pressure at the closed end can be 
expected to remain at an approximately constant low level as the detonation front proceeds through the 
chamber. The thrust is approximated during this phase by a linear relationship: 


F l=fo +at ’ 

where / 0 is a constant and a is the time rate of change in thrust. 


(60) 
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The timescale for phase I depends upon the blastwave characteristics and the engine geometry. In 
general, it is reasonable to assume that the thrust during phase II may be approximated by an exponentially 
decaying function: 


F u =Pex p 


~( r ~ T l) 

T II 


(61) 


where yS is a constant and T n is the characteristic decay time for blowdown. 

The model fits to the ensemble-averaged thrust profile are also shown in figure 25. It is interesting 
to estimate the average pressure at the closed end of the tube from the fitted value of/ 0 = 51 lb. In this case, 



511b 
0.785 in. 2 


~ 64 psi . 


(62) 


This implies a pressure ratio of only 4 with respect to the atmosphere. The actual average chamber pressure 
ratio is much higher of course when one averages over the entire centered expansion. 

The / for our apparatus was deduced from the ensemble-averaged thrust profile. This was done by 
first evaluating the total impulse of the system as 


/ = J'cyde F (t) dt = 0.43 lb -sec 


(63) 


The mass of propellant consumed per cycle m could be defined in terms of the scavenging ratio R s 


m p=l'™'’m(t)dl = V cPs R s . 

Assuming R s = 1, one deduces m = 0.982 g (2.16 x 10“ 3 lbm). The I sp is then determined by 


(64) 


sp W P (s!s c ) m p 

where w is the propellant weight, g is the Earth’s gravitational acceleration, and g c is a conversion factor. 
This yields an I sp of 200 sec. 

These thrust measurements alone are of great practical interest. Although the many potential 
advantages of PDE’s have been duly noted, the open scientific literature contains no meaningful comparison 
with conventional aerospace propulsion systems, and much confusion exists with respect to representative 
performance parameters. Specific impulse estimates have generally been inferred from comparisons of 
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steady Brayton and Humphrey cycles, but these estimates do not take into account the unsteady nature of 
the blowdown process and the dramatic manner in which the real process differs from the idealized 
thermodynamic cycle. 

In some cases, / comparisons have been based on the pressure rise created by the detonation 
process, but this is not the pressure that provides thrust in the PDE cycle. The centered expansion that 
follows a detonation in a closed tube reduces the detonation pressure ratio from =18 for stoichiometric 
H 2 /0 2 to an overpressure of =10, a value more representative of actual thrust performance. Measured / ’ s 
for PDRE’s have been reported in the range of 150 to 250 sec, but these have not been compared to 
conventional rocket engine performance at equivalent conditions. Furthermore, no I sp estimates have included 
the effects of a nozzle. 

From this viewpoint, it is useful to compare the thrust performance of the present pulse detonation 
apparatus with that of an equivalent conventional rocket engine. Equivalent conditions are defined as the 
same oxidizer-to-fuel ratio, the same nozzle type, and similar combustor pressure. Because the apparatus is 
based on a constant cross-sectional area with no divergent section, the conventional rocket performance 
should be computed based on a convergent nozzle only. Unfortunately, specification of chamber pressure is 
problematic due to the centered expansion process following the detonation wave. To provide perspective, 
a conventional rocket calculation was performed with the NASA SP-273 code using the C-J detonation 
pressure for C 2 H 2 /0 2 (33.6 atm) as the rocket chamber pressure. The computed I sp of 120 sec for the 
nozzleless rocket may be compared with the measured value of 200 sec for the pulse detonation apparatus. 
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4. CONCLUSIONS 


The prospects for realizing an integrated pulse detonation propulsion and MHD power system were 
examined. This included a critical evaluation of operational power requirements for direct detonation initiation 
of various fuel-oxidizer mixtures and the pumping power requirements for effective scavenging of the 
chamber. In addition, a series of experiments were conducted to investigate the basic engineering performance 
characteristics of a laboratory-scale, single-cycle pulse detonation-driven MHD generator. The experimental 
results reported here were only for stoichiometric mixtures of C 2 H 2 and 0 2 with an atomized spray of 
CsOH dissolved in alcohol as an ionization seed in the active MHD region. A plasma diagnostic channel 
was first attached to the detonation tube as a means of investigating ionization properties behind the detonation 
front. A continuous-electrode Faraday channel was subsequently attached to the detonation tube as a means 
of investigating MHD power generation characteristics. Time-resolved thrust characteristics of the single- 
cycle device were determined as well using a dynamic load cell incorporated into the thrust stand. 

The critical energy requirements for direct detonation initiation of various fuel-oxidizer mixtures 
were deduced from available experimental data and theoretical models. The resulting data were critiqued 
and summarized in tabular form. Using an established model, analyses were made of the variation in 
critical energy with igniter deposition power for stoichiometric H 2 /0 2 and H 2 /air mixtures. In general, the 
use of practical fuel-oxidizer mixtures in high firing rate PDE’s that may require a substantial amount of 
electrical power for direct detonation initiation was concluded. For H 2 /0 2 mixtures, the critical energy for 
direct detonation initiation ranged from 0.1 to 1 J/cm 2 , depending on the average deposition power; that is, 
the critical energy fell as the deposition power increased. The critical energy for H 2 /air mixtures was more 
than an order of magnitude higher than the H 2 /0 2 case. 

The pumping power required for effective scavenging and recharging of the chamber was also 
addressed. This was accomplished through the introduction of a scavenging ratio parameter and a scavenging 
efficiency parameter. These parameters were then used to develop scaling relationships for the required 
scavenging pressure ratio in the engine as well as the pumping power needed to achieve the required 
scavenging pressure. Representative calculations based on these scaling relationships demonstrate that 
significant pumping power may be necessary for effective scavenging of practical engines. 

In practice, this operational power must be obtained from either an independent, onboard power 
generation unit or from a power generator that is integrated into the engine. An independent power unit is 
feasible, but it negates the potential benefits of PDE’s associated with reduced weight, simplicity, and cost. 
On the other hand, PDE’s, unlike gas turbines, lack rotating shaft power to drive an integrated electrical 
generator of conventional design. If the electrical conductivity of the burned gas behind the detonation 
front is sufficiently high, however, it appears that an MHD generator may perform this function adequately. 

The laboratory-scale experiments of this work yielded data on detonation wave ionization properties, 
MHD power extraction characteristics, and thrust performance. Useful scaling relationships were also 
developed which permit extrapolation of data from small-scale units to practical-sized devices. 
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Measurements using a short-length plasma diagnostic channel attached to the end of the detonation 
tube confirmed the attainment of detonation conditions (p 2 IP\ ~ 34 and D ~ 2,400 m/sec) and yielded 
detailed information with respect to ionization quality and gas-dynamic structure. For instance, the measured 
effective electrical conductivity of the burned ionized gas was a ~ 6 S/m according to both two- and four- 
electrode probing techniques, and the measured dimensionless near-electrode voltage drop was A ~ 0.6. 
The conductivity value was low relative to theoretical expectations and indicated that the seeding mechanism 
of this experiment was far from optimized. The magnitude of the dimensionless near-electrode voltage 
drop is characteristic of small-scale devices and can be minimized by increasing electrode temperature and 
increasing channel height. The delayed appearance ( ~1.5 msec behind the detonation front) of a secondary 
conductive wave that was most likely associated with flow reversal in the channel during the latter portion 
of the blowdown process was also observed. This effect is due to the single-cycle operational limit of the 
device in which the head end of the tube remains closed throughout the cycle. It is also worth noting that a 
receding Cher-Kistiakowsky secondary compression wave emerged in the acquired static pressure waveforms 
at a time delay of =400 jusec behind the detonation front. 

Power extraction experiments were conducted using a short-length, continuous-electrode Faraday 
channel. The channel was placed in the bore of a magnet and attached to an active loading circuit in order 
to simulate higher effective magnetic induction. The following conclusions were reached with respect to 
power extraction: The measured peak open-circuit voltage was V oc = uBh - 10 V, yielding an effective 
burned gas velocity of 660 m/sec; measured load line characteristics were linear away from the extreme 
loading conditions and indicated decreasing internal plasma resistance with increasing magnetic induction; 
the energy density scaling parameter ranged from 10 to 10 3 J/m 3 as the effective magnetic induction was 
varied from 0.6 to 4.2 T; optimal power extraction shifted from a load impedance of 10 Q, when v 0 = 0 V, 
to 5 £2, when v 0 = 100 V, due to decreasing internal plasma resistance as the effective magnetic induction 
and Joule dissipation increased; theoretical predictions for the energy density scaling parameter were in 
reasonable agreement with measurements, given the uncertainties in parameter values; and by scaling up to 
a practical-sized device and limiting the near-electrode potential drop to 10 percent of the induced Faraday 
potential, a fivefold to tenfold increase in power generation performance with respect to our laboratory- 
scale results is anticipated. 

Time-resolved thrust measurements were also obtained and it was found that the thrust profile 
could be separated into two distinct phases. The first phase is characterized by a short-term, nearly constant 
thrust level while the detonation wave propagates through the chamber. The second phase is characterized 
by a decaying exponential-like thrust profile associated with chamber blowdown. The I sp of the single- 
cycle device was determined to be =200 sec. Calculations for an equivalent nozzleless rocket yielded an I 
of 120 sec. SP 
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